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We discuss the two-loop integrals necessary for evaluating massless 2 — » 2 scattering amplitudes. As a test 
process, we consider the leading colour two-loop contribution to qq — > q'q'. We show that for physical scattering 
processes the two Smirnov-Veretin planar box graphs I\ and I2 are accompanied by factors of 1/(D — 4) thereby 
necessitating a knowledge of both Ii and I2 to 0(e). Using an alternative basis I\ and the irreducible numerator 
integral 73, the factors of 1/(D — 4) disappear. 



1. Introduction 

Two-to-two scattering processes are well known 
to be one of the most basic probes of the funda- 
mental interactions of nature. In hadron-hadron 
collisions, parton-parton scattering to form a 
large transverse momentum jet tests the point- 
like nature of the partons down to distance scales 
of 1CU 17 m. However, extracting useful results 
from experimental data requires both plentiful 
data and accurate theoretical calculations. For 
example, the single jet inclusive transverse en- 
ergy distribution observed by the CDF collabora- 
tion in Run I at the TEVATRON indicated pos- 
sible new physics at large transverse energy [ [j]. 
Data obtained by the DO collaboration [ || was 
more consistent with theoretical next-to-leading 
order expectations, however, because of both the- 
oretical and experimental uncertainties no defi- 
nite conclusion could be drawn. The experimen- 
tal situation may be clarified in the forthcoming 
high statistic Run II starting in 2001. The the- 
oretical prediction may be improved by includ- 
ing the next-to-next-to-leading order perturba- 
tive predictions. This has the effect of (a) reduc- 
ing the renormalisation scale dependence and (b) 
improving the matching of the parton level theo- 
retical jet algorithm with the hadron level experi- 
mental jet algorithm because the jet structure can 
be modelled by the presence of a third parton. 
Varying the renormalisation scale up and down 
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by a factor of two about the jet transverse energy 
leads to a 20% (10%) renormalisation scale un- 
certainty at leading order (next-to-leading order) 
for jets with Et ~ 100 GeV. The improvement in 
accuracy expected at next-to-next-to-leading or- 
der can be estimated using the renormalisation 
group equations together with the known leading 
and next-to-leading order coefficients and is at 
the 1-2% level. Of course, the full next-to-next- 
to-leading order prediction requires a knowledge 
of the two-loop 2^2 matrix elements as well as 
the contributions from the one-loop 2 — > 3 and 
tree-level 2^4 processes. 

In this talk, we wish to review the re- 
cent progress that has been made towards the 
analaytic evaluation of the two-loop matrix ele- 
ments relevant for massless 2 — > 2 scattering. As 
can be seen from Table 1, the number of Feyn- 
man diagrams contributing to the basic parton 
scattering processes increases dramatically with 
the number of loops. The one-loop graphs are 
those computed by Ellis and Sexton [ f| in 1986. 
The much more numerous two-loop graphs may 
be either products of one- loop graphs, self-energy 
insertions or genuinely new topologies. It is the 
latter class which has proved to be a major stum- 
bling block. However, in the last twelve months, 
all of the necessary integrals have been computed 
and a complete basis set of master integrals now 
exists. We note in passing that a particular two- 
loop helicity amplitude for gg — > gg scattering 
has been calculated by Bern, Dixon and Kosower 
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Process 


Tree 


One loop 


Two loops 


99 55 


4 


72 


1531 


99 -> 99 


3 


29 


563 


qq -> q'q' 


1 


10 


186 



Table 1. Numbers of Feynman diagrams con- 
tributing to 2 — > 2 parton scattering processes 



2. Master Integrals 

The complete set of massless master integrals 
comprises the trivial topologies of single scale 
integrals which can be written as products of 
Gamma functions, 

-e-w -oo-w -cn> 

the less trivial non-planar triangle graph [ B| , 



(s) 




two scale integrals that are related to the one-loop 
box graphs [ fm , 



the planar double boxes [ 0, ^] 
(«>*) 



(s,t) 



(s,t) 



which we denote I\ and I2 respectively together 
with the non-planar double boxes [ ^, [l(| , 



The Mandelstam variables s and £ represent the 
kinematic scales involved in the integral while 
the blobs on the propagators represent an ad- 
ditional power of that propagator. The latter 
blobbed graphs are necessary to evaluate tensor 
integrals. In other words, starting from a planar 



or non-planar box tensor integral, it is not pos- 
sible to reduce the powers of all propagators to 
unity and the second master integral is required 
[ ||, [l0|. The scalar planar [ fjj and non-planar [ 
^ integrals themselves were evaluated as multiple 
Mellin-Barnes integrals and represent significant 
achievements in the field of Feynman diagram- 
mology. 

3. Application: qq — > q'q' 

As a precursor to a full two loop calculation 
of all massless 2 — > 2 scattering matrix ele- 
ments let us concentrate on the specific process 
q(pi)q(p2) — * Q 1 (P3)q' {pi) where the lightlikc mo- 
mentum assignments are in parentheses. 

The amplitude M has the perturbative expan- 
sion, 

M = g 2 M + gtMi + g 6 s M 2 + 0(g s s ), 

in terms of the tree level (Mo), one- loop (Mi) 
and two- loop (M2) amplitudes. The squared and 
summed tree-level amplitude is given by, 



\M \ 2 = 2(N 2 



(i) 



where s — (P1+P2) 2 , t = (P2—P3) 2 and u = —s—t. 
The one- loop amplitude Mi was first calculated 
by Ellis and Sexton [ ||| and contributes to the 
cross section at O(gf), or next-to-leading order. 
The two- loop amplitude M2 first contributes at 
0(g®) through its interference with the tree level 
amplitude Mo and has the following colour struc- 
ture, 



M 2 M\ +M\M = (N 2 - 1) 
1 



AN +B + C 



N 2 



DN F N + E- 



,Nf 
~N 



where N is the number of colours and Nf the 
number of light fcrmions. The leading colour am- 
plitude A is gauge invariant and contains only 
planar diagrams such as those shown in Fig- 
ure 1. The amplitudes suppressed by powers of N 
(B, C, D and E) contain the non-planar graphs. 
As a first step in carrying out the analytic eval- 
uation of the two loop graphs, we therefore focus 
on the leading colour amplitude A. 
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Figure 1: The planar box graphs for qq 
contributing at leading colour 



q'q' 



4. Auxiliary diagram 

To handle all possible permutations of planar 
diagrams it is convenient to work with the auxil- 
iary diagram shown in Figure 2, 



Pi 



Pi 



P2 




P2 



Figure 2: The general auxiliary planar diagram 
with propagator i raised to the power vi 
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where A x = fcf, A 2 = (fci+pi) 2 , A 3 = (fci+pi 2 ) 2 , 
A 4 = (fci+pi 23 ) 2 , A b = kl A e = (k 2 +pi) 2 , A? = 
(k 2 +Pu) 2 , A s = (fe 2 +P123) 2 , A 9 = (fc 2 - fci) 2 
and p y = ]3j + pj and p ljk = Pi + Pj + Pk- The 
ith propagator is raised to the power i>i. Scalar 
integrals have all Vi = 1 or 0. For example, 



(s,t) = 7^(1,1,1,0,1,0,1,1,1) 
(s,t) = 7^(1,1,0,0,0,0,1,1,1) 
Q (s,t) = I D (1, 1,1, 0,0, 0,0, 1,1) 



By permuting the arguments, we obtain the other 
orientations. For example, 



(t,s) = 7^(1,1,0,1,0,1,1,1,1). 



For the interference of tree-level with two-loop 
graphs, loop momenta in the numerator are al- 
ways contracted with either external or loop mo- 
menta. These dot-products can always be written 
as combinations of the 9 propagators so that ten- 
sor integrals appear as generalised scalar integals. 
Negative values of Vi correspond to irreducible nu- 
merators. For example, the planar box integral 
with one irreducible numerator on the left hand 
loop can be written, 



® 



= 7°(1,1, 1,-1, 1,0, 1,1,1) 



5. General procedure 

The general procedure for computing the two 
loop graphs is as follows: 

1. Use QGRAF [ |ll| to generate the Feynman 
diagrams 

2. Multiply by tree- level and compute traces 

3. Identify combinations of scalar and tensor 
auxiliary integrals 

4. Exchange tensor integrals (yt < 0) for aux- 
iliary integrals in higher dimension with 
higher powers of propagators (z/, > 0) [ [l2| 

5. Apply integration-by-parts (IBP) identities 
[ [l3| to reduce general auxiliary integrals 
to combinations of the Master Topologies 

6. Apply specific IBP identities to reduce Mas- 
ter Topologies to Master Integrals (vi = 1) 

6. Planar box graphs 

However, because the specific IBP identities 
for the planar box are quite complicated [ || , we 
choose an alternative method for the tensor pla- 
nar box graphs. In particular we try to stay in 
D ~ 4. To do this we adopt the approach of 
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Gehrmann and Remiddi [ [l4|] . The idea is very 
simple both in concept and in implementation. 
We characterise a loop integral with three num- 
bers, t the number of different propagators in the 
denominator, r the sum of powers of propagators 
in the denominator and s the sum of powers of 
propagators in the numerator. For example, 



© 



(M) 
(s,t) 
(s,t) 



t = r = 7, s = 



7, 



t = 7, 







When acting on a loop integral It, r ,s, the IBP 
(and Lorentz invariance [ [l4j) identities produce 
more complicated integrals with the same topol- 
ogy I t} r+i, s an d It,r+i,s+i- simpler integrals with 
the same topology It. r -i, s and It, r -i,s-i as well 
as simpler topologies It-i, r ,s- By applying each 
identity to each It,r,s , we can form a linear system 
of equations from which the more complicated in- 
tegrals can be eliminated. 

For example, when t — 7, the number of in- 
tegrals for a given value of r and s is shown in 
Table 2. The two master integrals [ |[ || have 





s 







1 


2 


3 


4 




7 


1 


2 


3 


4 


5 


r 


8 


7 


14 


21 


28 


35 




9 


28 


56 


84 


112 


140 




10 


84 


168 


252 


336 


420 



Table 2: Numbers of integrals with different val- 
ues of r and s for t = 7 (taken from [ [T3] ) . 



t = 7 and r = 7 and r = 8 respectively. Tensor 
integrals (corresponding to r — 7 and s > 0) lie 
on the first row of Table 2. Using only these inte- 
grals as seeds for the IBP identities and eliminat- 



ing the unknowns using linear algebra, we imme- 
diately obtain all tensor integrals for the planar 
box in D ~ 4 in terms of 1\ and the irreducible 
numerator graph 



(rather than 1\ and I2) together with simpler 
pinched integrals that can be straightforwardly 
simplified. If we denote the s = i + j,t = r = 7 
planar box integral as 



I D (1, 1, -J, 1,1,1) 



then, 
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and, for example, the second rank tensor integrals 
are given by, 



©© 



st 
~2 



M) 



'is 
~2 



© 



8{D - 3) 

(£>-4) 
(7s + 9t) 



(s,t) 

7Qm) 

17(£>-3)(3Z>-10) _0~ (t 

-e-e 



+ 



2s(D-4) 2 
2(3£> - 8) (3D - 10) 

s 2 (L>-4) 3 
9(3£>-10)(3D-8)(£>-3) 
+ st{D - 4) 3 

2{D - 3)((2L> - 5)s + 2(D - 3)t) 
s 2 {D~Af 

_ {D-A)st 
~ 2{D - 3) 



-e-w 

-00 w 



(3(D - 4)s - 2i) 
2{D - 3) 

Q~(«,t) 
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+ 



+ 



(7s + &)(£> -4) I A , 

s{D-3) Y\ M) 

(13 3 -2*)(3P-10) _TTs- , , 
2^4) CD_( fi ) 

2(s + 2i)(3D-8)(3D- 10) 

2s 3 (D-3)(D-4) 2 

9(3D - 8) (3D - 10) 



st(D - 4) 2 ©~ (<) 

2((2D-7) S + 2(D-4)t) r^r^ ( \ 

We see that the coefficients of both I\ and I3 are 
always finite as D — ► 4 and that the IBP identities 
have not introduced fake singularities. 

7. Leading Colour Matrix Elements 

In D dimensions, the leading colour two-loop 
amplitude A for qq — * g'g' scattering has the fol- 
lowing structure, 

9 / 

A = - 



(4tt) j 



-«3 



M) 



(*.*) 



«2 
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(«>*) 



«4 © 



(t, a) 



no 



He- 



+a 5 

+a 7 \/\ (s.t) 
+on Q («) 

+«i 3 -00~ (s) 



We have calculated the coefficients ai — 014 in 
arbitrary dimension. For example, the first two 
coefficients of the planar box graphs are given by, 




01 = 



(llD 2 -53D + 54) 

4(D - 3) 
(J 2 + 2D - 18) 2 2 
(D-3) " 



s J £ - 4 sf 3 



02 = 



(15D 2 - 71D + 69) 
2(D-3) 



2(7D - 24) 2 
(D-3) 



(D-3) 

We see that both a% and 02 are well behaved as 
D — > 4 indicating again that the IBP identities 
have not introduced fake singularities. 

8. Relationship between Master Integrals 

We can also use the Gehrmann-Remiddi ap- 
proach [ [lJJ to find a relation between the irre- 
ducible numerator master integral I3 and those 
of Smirnov and Veretin (Ii and I 2 ) [ @, Sup- 
pressing the simpler pinched integrals, we find, 



© 



(3D - 14)s 
2(D-4) 



(D - 6)st 



2(D-4)(D-5) 
pinchings. 



The presence of the D — 4 factor in the denomi- 
nator immediately indicates a problem. Both I\ 
and I2 have 1/e 4 leading poles and it would ap- 
pear that I3 ~ 1/e 5 . This is not the case as close 
examination of I\ and I2 shows that in this com- 
bination the 1/e 4 poles and descendents cancel 
completely so that 
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as we expect. However, the finite parts of I3 are 
controlled by the 0(e) parts of I\ and 12- 

I3 can also be written in terms of derivatives of 

h, 



© 



(s(D-5)-t) 



(D 
id 



4) 

d_- 

dt. 



(D-4) 
+ pinchings 

which again indicates that the 0(e) part of I\ is 
necessary to determine the 0(e°) part of I^- In 
fact this additional part is not very difficult to 
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obtain either by considering the differential equa- 
tions for Ii and I 2 at t = — s [ [l5| or by explicit 
evaluation of the Mellin-Barnes integrals [ |l6| . As 
expected, the e expansion for ^3 through to O(e ) 
contains quadrilogarithms at worst. 

9. Results 

Using the analytic expansions around D ~ 4 
we can evaluate the leading colour two-loop am- 
plitudes for qq — > q'q'. Expanding both the in- 
tegrals and coefficients a\ - an we find that the 
leading singularities are proportional to tree level, 



m 2 mI+m\Mq = 

AN 2 L(l + e) 2 L(l-e) 4 1 



(4ir) D 



r(l-2e) 



-2c 



The leading pole has the same coefficient as the 
square of the one-loop amplitude. The remaining 
pole and finite contributions are in the process of 
being checked. 

10. Summary 

We have made a study of two-loop amplitudes 
of massless 2-^2 scattering by considering qq — > 
q'q' as a trial process. At leading colour, only pla- 
nar graphs contribute and we have expressed the 
amplitude as a sum over the basis set of two-loop 
master integrals. It turns out that for this pro- 
cess, in reducing the tensor integrals to scalars, 
factors of 1/(D — 4) are generated multiplying 
both of the Smirnov-Veretin planar box graphs, 
l\ and These factors do not cancel in the 
physical process thereby necessitating a knowl- 
edge of both 1 1 and I2 to 0(e). Alternatively, if 
one uses the basis 1\ and the irreducible numera- 
tor integral I3 , the factors of 1/ (D — 4) disappear. 
Of course, evaluating 23 also requires the 0(e) 
of 1\ but this has now been calculated [ [l5L jig . 
We therefore have the ingredients to evaluate the 
leading colour two-loop amplitude for qq — > q'q' 
for D ~ 4 and the leading singularities agree with 
expectations. 
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